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Abstract 

We are carrying out a series of photometric monitoring to measure the rotation 
periods of members in the young (3 Pictoris Association, as part of the RACE- 
OC project (Rotation and Activity Evolution in Open Clusters). In this paper, we 
present the results for HD 155555C which is believed to be physically associated 
to the spectroscopic binary V824 Ara (HD155555) and thus constituting a triple 
system. We collected B, V, and R-band photometric data timeseries and discovered 
from periodogram analysis the rotation period P = 4.43 d. Combined with stellar 
radius and projected rotational velocity, we find this star almost equator-on with an 
inclination i ~ 90°. The rotational properties of HD155555C fit well into the period 
distribution of other (3 Pic members, giving further support to the suggested mem¬ 
bership to the association and to its physical association to V824 Ara. A comparison 
with Pre-Main-Sequence isochrones from various models allows us to estimate an 
age of 20±15 Myr for this triple system. 

Key words: Stars: activity - Stars: low-mass - Stars: rotation - Stars: starspots - 
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1 Introduction 


We are carrying out a systematic investigation aimed at determining the rota¬ 
tion periods and the photospheric activity levels of all confirmed and candidate 
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members of the young /3 Pictoris stellar Association (Zuckerman et al. 2001). 
This investigation is part of the RACE-OC project (Rotation and Activity 
Evolution in Open Clusters; Messina [ 2007] ). In an earlier study by Messina 
et al. pOTO] . we focused on only conhrmed members and derived the rota¬ 
tion periods for a total of 27 of them either from the literature or from our 
analysis of photometric timeseries in the ASAS (All Sky Automated Survey; 
Pojmanski [1997] ) and Super WASP (Wide Angle Search for Planets; Butters 
et al. [ 2010] ) archives. However, a few members have remained with unknown 
or uncertain periods for a number of reasons: they were not in the archives; 
the low-photometric precision in combination with low variability level pre¬ 
vented us from rotation period detection; they were in spatially unresolved 
multiple systems. Meanwhile, new candidate members have been proposed, 
e.g., by Lepine et al. [ 200^] . Kiss et al. [ 2011] . Shkolnik et al. [2012 ]. Schlieder 
et al. [2UT2] . Malo et al. [2014] . In order to have a complete rotational charac¬ 
terization of the association, we have undertaken a number of collaborations. 
In the present paper, we present the results of the monitoring project carried 
out at the York Creek Observatory on HD 155555C, one of the first proposed 
members to this Association. 

We focus our investigation on the pre-main-sequence triple system HD 155555 
consisting of a SB2 spectroscopic binary, HD 155555AB, and a nearby single 
dM component, HD 155555C. HD 155555AB has been identified as member 
of the young stellar Association f3 Pictoris by Zuckerman et al. [2001] , The 
nearby dM star has been suggested to be physically connected to the SB2 
by Pasquini et al. [ 1989] and Pasquini et al. [1991 ]. and therefore to be a 
member of the same Association. However, no parallax of the latter is mea¬ 
sured, and membership to this stellar system as well as to the Association are 
only inferred, though significantly supported by a number of evidences. Our 
rotational investigation aims to add additional information to get a better 
characterization of this dM component. 

For instance, this system, like other multiple systems, see, e.g., Messina et al. 
[20 Hj, is particularly interesting since a few parameters are equal, like age, 
initial chemical composition, and mass for some systems. These commonali¬ 
ties allow us to explore the dependence of their rotational evolution on other 
parameters, such as initial angular velocity and environment properties (pres¬ 
ence of disc or nearby stcllar/planetary companions), shedding light on the 
causes of dispersion observed in the rotation period distributions of all young 
clusters/associations. 


2 Literature information 


HD 155555 is a triple stellar system in the southern constellation of Ara at 
a distance d = 31.4±0.5 pc (van Leeuwen [20U7 ]). This system consists of a 
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double-line spectroscopic binary HD 155555AB and of a fainter M dwarf com¬ 
panion HD 155555C at an angular distance of about 33", which corresponds 
to about 1040 AU. 

The spectroscopic binary HD 155555AB was discovered by Bennet et al. [ 1967] 
to be composed of G5IV + KOIV stars that rotate with an orbital period P 
= 1.687 d. First evidence of photometric variability of this SB2 was found by 
Eggen [ 1978] after comparing his photometry collected in 1978 (V = 6.83, 
B—V = 0.835, U—B = 0.29, R = 6.46, and R—I = 0.325 mag) with the earlier 
photometry collected by Stoy [1963] (V = 6.67 and B—V = 0.80 mag), even 
though the latter also included the dM companion, owing to the large aper¬ 
ture of 44" used to collect the photoelectric photometry. Subsequent long-term 
photometric monitoring by Cutispoto (see, e.g., Cutispoto 1998 and references 
therein) revealed robust evidence of the existence of short- and long-term pho¬ 
tometric variability. HD 155555AB has the variable’s name V824 Ara in the 
Combined General Catalogue of Variable Stars |Kholopov 199 81. The most 
recent determination of orbital and physical properties, and first magnetic 
maps of both components were provided by the spectro-polarimetric study 
of Dunstone et al. [ 1998] , Most recent surface temperature maps have been 
obtained with an improved Doppler imaging technique by Kristovics et al. 
[ 2013] . The membership of HD155555AB to the f3 Pictoris moving group was 
first suggested by Zuckerman et al. [2001 ] based on distance, UVW velocity 
components, high usini, and Lx/ Lboi ratio. Also the photometric variability 
observed by Eggen (1978) and Cutispoto (1998) are in agreement with the 
youth of the system. 

The present investigation is focused on the fainter HD 155555C component 
(also named LDS 587B) which is a M3.5 dwarf (Pasquini et al. 1991). Eggen 
[ 1978] reported the following magnitudes and colors: V = 12.82, B—V = 1.54, 
U—B = 1.05, R = 11.40, and R—I = 1.205 mag, noticing both UV and IR flux 
excess with respect to standard colors of main sequence stars. This component 
is a strong X-ray source (EXO 171224-6653.9), which was a serendipitous dis¬ 
covery by Barstow et al. [ 1987] during EXOSAT observations of HD 155555AB 
(EXO 171217-6653.8). Pasquini et al. [1989] identified this X-ray source with 
HD 155555C, and subsequent spectroscopic investigation (Pasquini et al. 1991) 
showed high levels of magnetic activity, the Hel, Nal D2 and D1 lines all be¬ 
ing in emission. The Ha line was also found to be in emission and variable 
from night to night with the equivalent widths ranging from 6.2 to 9.7 A. 
Pasquini et al. [ 1991 ] could measure only an upper value (< 3 kms -1 ) of the 
projected rotational velocity. Their study concluded that HD155555AB and 
HD155555 C are physically associated and are likely PMS stars. They noticed 
that radial velocities are similar (V r = 2.74kms -1 for HD 155555AB com¬ 
pared to V r = 3.5kms -1 for HD 155555C. Torres et al. [ 2006] also find that 
the UVW velocity components and radial velocity of HD155555C are similar 
to those of HD155555AB, suggesting their physical association. 

The physical association is also supported by proper motion studies by Mar¬ 
tin et al. [ 1995] . who find the angular separation changed from 33" (about 
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Fig. 1. Light curves of HD155555C phased with the rotation period in three filters. 

Blue bullets denote observations collected after HJD 2456761 to better show the 
light curve evolution. The rectangle denotes a possible flare event. 

in 1929) to 32.6" in 1994, suggesting similar proper motions. Similarly, the 
Washington Visual Double Star Catalog [Mason et ah 200T] reports a small 
variation from 33" in 1920 to 34" in 2000. In the literature, two other measure¬ 
ments of projected rotational velocity are available from Torres et al. | 2006 j 
usini = 6.0 ±1.2 km s^ 1 and from Weise et al. 1 2010 r;sinf = 8.8 ±0.9 kms^ 1 . 

This system has been searched for possible planetary companions using the 
NACO adaptive optics system at the Very Large Telescope (VLT) [Masciadri et ah 2005j . 
and the MMT telescope [Biller et al. 2007j . It was included in the SIM Plan- 
etQuest Key Project [Tanner et al. 2007] , and more recently in the Gem- 
ini/NICI planet-finding campaign [Biller et al. 2013] . None of these investi¬ 
gations have so far detected the presence of any additional companion. 


3 Photometric observations 


HD 155555C was observed in 2014, from Feb 28 to May 10, for a total of 21 
nights. Observations were collected at the York Creek Observatory (41° 06' 06"S; 
146° 50' 33"E, Georgetown, Tasmania) using a f/10 25cm Takahashi Mewlon 
reflector, equipped with a QSI 683ws-8 camera, and BVR standard Johnson- 
Cousins filters. The telescope has a FoV of 24.5'x 18.5', and the plate scale is 
0.44"/pixels. On each night observations were collected from 3 to 6 consecu- 
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Fig. 2. Correlations among magnitude variations in the three filters. Solid lines are 
weighted linear fits, whereas labels indicates the linear correlation coefficient (r) and 
the fit slope with its uncertainty. 
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Fig. 3. left panel : Lomb-Scargle periodogram of HD155555C (solid line) with indica¬ 
tion of the power peak corresponding to the rotation period P = 4.43 d. The dotted 
red line indicates the window spectral function, whereas the horizontal dashed line 
the power level corresponding to a FAP = 0.01. right panel: Clean periodogram. 

five hours, using 120-s integrations in each filter. A total of 360 frames in B, 
359 in V, and 400 in R filter were obtained (see Table]!]). 


After bias subtraction and flat fielding, we extracted the magnitude timeseries 
for a total of 30 stars. These series were cleaned by applying a 3a threshold 
to remove outliers. The differential photometry was performed using the tech¬ 
nique of ensemble photometry (Gilliland et al. 1988; Bailer-Jones et al. 2001; 
Everett et ah 2001). The differential magnitude of the target star is computed 
with respect to the average magnitude of a large number of non-variable ref- 
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Fig. 4. Spectral Energy Distribution of HD155555C. 

erence stars. In this averaging process, the uncertainties of the ensemble star 
magnitudes due to statistical fluctuations as well as short-term small inco¬ 
herent variations will cancel each other. The uncertainty in the magnitude 
of the artificial comparison will therefore be smaller than the uncertainty on 
the magnitude of a single star and this, in turn, will produce less noisy light 
curves. 

We started with 30 sufficiently bright stars, which were isolated, distributed 
all over the frame but not found to be close to the CCD edges, and common to 
all the frames. Having selected the ensemble, we computed (i) the magnitude 
of the artificial comparison star by averaging the instrumental magnitudes 
of ensemble stars, (ii) time-series differential magnitudes for each star of the 
ensemble and (iii) mean, median and the standard deviation (a) associated 
with each time series. While constructing the final ensemble, we excluded those 
ensemble stars whose standard deviation was larger than the threshold a value 
that we fixed to 0.010 mag. After a number of iterations, we were left with 
six stars whose standard deviations were smaller than 0.010 mag. In Table[2] 
we list the 2MASS designation jCutr i et ah 2003] comparison stars, the nearby 
optical sources from USNO catalogue [ Monet et al. 2013 ]. the coordinates and 
the J-band magnitude. These stars were used to build the ensemble comparison 
star and to obtain the differential magnitude time series of HD 155555C in B, 
V, and R filters. 

In Fig.[lJ we plot the B, V, and R light curves phased with the P = 4.43 d 


6 

















Table 1 

Log of observations. 
Observation # frames 


Date 

B 

V 

R 

28/02/2014 

8 

10 

8 

02/03/2014 

10 

10 

43 

19/03/2014 

4 

4 

34 

24/03/2014 

3 

3 

3 

25/03/2014 

19 

16 

17 

03/04/2014 

3 

3 

3 

04/04/2014 

5 

6 

5 

05/04/2014 

20 

24 

24 

06/04/2014 

22 

20 

14 

12/04/2014 

22 

22 

19 

14/04/2014 

13 

15 

15 

24/04/2014 

32 

32 

29 

27/04/2014 

17 

17 

18 

30/04/2014 

33 

33 

33 

03/05/2014 

24 

26 

26 

04/05/2014 

20 

18 

17 

05/05/2014 

18 

18 

17 

06/05/2014 

27 

26 

27 

07/05/2014 

25 

24 

19 

08/05/2014 

10 

10 

9 

10/05/2014 

25 

22 

20 


rotation period (see Sect. 5). The average photometric uncertainties associated 
to each magnitude are 0.014, 0.007, and 0.005 mag in the B, V, and R filters, 
respectively. 


4 Analysis 


As shown in Fig.[2j the photometric timeseries show that HD 155555C varies 
in all B, V, and R bands. The magnitude variations are strongly correlated to 
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Table 2 

Comparison stars. 


2MASS 

USNO-Bl.O 

RA 

DEC 

J mag 

2MASS J17174644-6656455 

09064-02437 

17 17 46.44 

-66 56 45.60 

8.696 

2MASS J17174439-6657287 

0230-0693233 

17 17 44.37 

-66 57 28.84 

11.285 

2MASS J17175417-6657550 

0230-0693339 

17 17 54.17 

-66 57 55.06 

12.085 

2MASS J17172719-6700468 

0229-0763917 

17 17 27.19 

-67 00 46.86 

11.457 

2MASS J17172092-6701280 

0229-0763856 

17 17 20.93 

-67 01 28.08 

9.411 

2MASS J17170984-6700057 

0229-0763745 

17 17 09.84 

-67 00 05.77 

12.383 


each other. The Pearson linear correlation coefficients are all positive, rev = 
0.85, r VR = 0.92, and r BR = 0.89, and all with a significance level > 99.9%. 
Weighted linear fits to the magnitude variations (solid lines in Fig. [2]) give the 
following slopes AB/AV = 0.76±0.11, AV/AR = 1.12±0.06, and AB/AR 
= 1.23±0.17. This is the behavior typically exhibited, according to Messina 
[2008]. by color-correlated and reddening stars, whose color variations arise 
from cool spots. Such temperature inhomogeneities generally give rise to light 
curves whose amplitudes (peak-to-peak) increase towards shorter wavelengths, 
that is the B light curve has amplitude larger than the V light curve that is 
larger than the R light curve. In our specific case, we note that whereas V 
(whose amplitude is AV = 0.07 mag) and R (whose amplitude is AR = 0.05 
mag) curves follow this trend, the B-band light curve has amplitude com¬ 
parable, if not slightly smaller (AB = 0.06 mag), than V (AB/AV < 1). 
Unfortunately, our B-band observations could not cover the rotation phases 
around the lightcurve maximum (phases 0.15-0.35). Therefore, the true am¬ 
plitude of variation in the B band may be larger than measured and certainly 
further observations will address this aspect. The V and R light curves, for 
which we have the longest timeseries, exhibit evidence of a significant varia¬ 
tion, the amplitude has decreased and the phase of light maximum migrated 
by A0 = 0.2, as likely consequence of the active regions growth and decay 
(ARGD), which are responsible for the observed variability. 


5 Rotation period search 


We used the Lomb-Scargle |Scargle 1982| and Clean [Roberts et al. 1987] pe- 
riodogram analyses to search for significant periodicities in the HD155555C 
magnitude timeseries related to stellar rotation period. The results of our anal¬ 
ysis are plotted in Fig.[3j In the left panel we plot the normalized Lomb-Scargle 
periodogram (solid line) with overplotted the spectral window function related 
to the data sampling (red dotted line). The horizontal dashed line indicates 











the power level corresponding to a False Alarm Probability (FAP) of 1%, that 
is the probability that a power peak of that height simply arises from Gaussian 
noise in the data. The FAP was estimated using a Monte-Carlo method, i.e., 
by generating 1000 artificial light curves obtained from the real one, keeping 
the date but scrambling the magnitude values. In the right panel we plot the 
Clean periodogram where the power peak arising from the light rotational 
modulation dominates, whereas all secondary peaks, arising from the aliasing, 
are effectively removed. Our analysis finds that the stellar rotation period is P 
= 4.43±0.05d. The uncertainty of the rotation period is computed following 
the prescription of Lamm et al. 2001 . 


6 SED analysis 


We used the available optical, near-IR, and IR photometry to build the ob¬ 
served spectral energy distribution (SED). The UBVRI magnitudes from Eggen 
|1978j . JHK magnitudes from 2MASS |Cutri et al. 2003] . and W1-W4 magni¬ 
tudes from WISE |Cutri et al. 2013j are listed in Table[3j The SED was fitted 
with a grid of theoretical spectra from the NextGen Model [Allard et al. 2012] 
and the best fit is obtained with a model of T e ff = 3300±100 K, log g = 4.0±0.3 
dex and metallicity [Fe/H] = 0.0 (see Fig|4j) . We note that no evidence of IR 
excess was found. 


7 Target parameters 


We can use effective temperature and luminosity to derive the stellar radius 
and the inclination of the rotation axis. We assume that HD155555AB and 
HD155555C are bound and, therefore, they have the same distance d = 31.4 
pc as measured by Hipparcos for HD155555AB (Zuckerman et al. 2004). Using 
the brightest observed magnitude V = 12.71 mag (Zuckerman et al. 2001), and 
the bolometric correction BCy = —2.03 mag from Pecaut & Mamajek [ 2013] . 
we infer the bolometric magnitude Mboi = 8.19±0.05 mag, the luminosity L = 
0.042±0.005 L 0 , and the stellar radius R = 0.65±0.05 R 0 . From these values 
and the stellar rotation period we infer the inclination of the stellar rotation 
axis i ~ 90° (adopting < usini > = 7.6 kms -1 ). 

In Fig. [5} we compare T ( . fi - and luminosity of all components of HD 155555 with 
a set of isochronones taken from Baraffe et al. |1998j . Siess et al. [ 2000j . and 
D’Antona & Mazzitclli [ 1997] , for solar metallicity to infer the age of each com¬ 
ponent and see if their coevalness hypothesis is sustainable. The physical pa¬ 
rameters of HD 155555A and B are taken from Strassmeier & Rice 2UUT) . with 
the only difference that luminosities are computed using the BC from Pecaut 
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Table 3 

Photometry from the literature. 


u 

B 

V R 

I 

J 

H 

K 

W1 

W2 

W3 

W4 

(mag) 

(mag) 

(mag) (mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

15.41 

14.36 

12.82 11.40 6 

10.19 

8.542 

7.916 

7.629 

7.528 

7.358 

7.204 

6.940 


& Mamajek [ 2013 ]. We see that the D’Antona & Mazzitelli [1997] model pro¬ 
duces the smallest age dispersion among all components with an estimated 
age of 15±7 Myr for the system and a mass M=0.25 M© for HD 155555C. The 
Baraffe et al. |1998j and Siess et al. [ 2000 ] models produce results that empha¬ 
size a larger dispersion between the age of the SB2 components and the age of 
HD155555C. According to these models HD155555AB has an age Aearaffe = 
30±10 Myr and Asi eS s = 24±8 My that are about a factor 2 older than the age 
of Asaraffe = 13±7 Myr and Asi ess = 12±4 Myr of HD155555C. Despite these 
systematic differences between the Siess et ah [ 2000] and Baraffe et ah [1998 ] 
models and the D’Antona & Mazzitelli [ 1997] model, the uncertainties on the 
age are quite large and our results are consistent with the age range reported 
by Zuckermann et ah [ 2001] who found the system components to be fitted 
by a set of isochrones in the age range from 5 to 30 Myr. Thus, the coevalness 
hypothesis is still sustainable with an age of 20±15 Myr. Moreover, such an 
age is also consistent with the age of 21±9 Myr that Mentuch et ah [2008 ] 
estimate for the f3 Pic Association and based on the Li depletion. We note 
that the low Li content (EWli < 20 mA) measured by Martin et ah [1995 ] in 
the component C is consistent with the values measured on other Association 
members of similar spectral type. On the contrary, the hotter components A 
and B are in the un-depleted part of the Li distribution, which prevents us 
from using the Li content to constrain their ages. 


8 Spot model 


To show that the presence of cool spots is the most typical explanation of the 
light modulation observed in a low-mass star as HD155555C we performed a 
spot modeling of the light curves. We modeled the observed multi-band light 
curves using Binary Maker V 3.0 [Bradstreet &; Steelman 2004] to infer some 
information on the geometrical and physical properties of spots on the pho¬ 
tosphere of HD 155555C. Binary Maker V 3.0 models are almost identical to 
those generated by Wilson-Devinney program Wilson & Devinney 1971 and 
uses Roche equipotentials to create star surfaces. In our modeling the second 
component is essentially “turned off” (i.e., assigned a near zero mass and lumi¬ 
nosity) in order to model a single rotating star. The gravity-darkening coeffi¬ 
cient has been assumed v = 0.25 jKopal 1959 , and limb-darkening coefficients 
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evolutionary track for a mass M=0.25 Mg. Models are from Baraffe et al. [1998] 
(top panel), Siess et al. [ 2000j (middle panel), and D’Antona & Mazzitelli [ 1997 ] 
(bottom panel). 

from Claret et al. [2012 ] were also adopted. We adopted T e g = 3300 K, and 
inclination i = 90° as input parameters. To find a suitable spot configuration 
we took into account that the light curve shape is asymmetric, the decreasing 
flux branch is shorter than the increasing flux branch. This suggests the exis¬ 
tence of two major spot groups separated in longitude. Unfortunately, we had 
no observations to be fit around the light curve maximum. 

Thanks to the dependence of the light curve amplitude on the photometric 
band wavelength, we could constrain the spot temperature contrast and better 
determine the area of the spots responsible for the flux rotational modulation. 
The temperature contrast between spots and surrounding photosphere was 
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Fig. 7. Rotation period distribution of /3 Pic members from Messina et al. mm- 
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12 

















found T sp0 i/T p hot = 0.85, which is a value generally measured in dM3 spot¬ 
ted stars (see, e.g., Berdyugina 2005). The only free parameters in our model 
remained the spot areas and the spot longitudes. The spots latitude cannot 
be constrained by photometry alone, especially in the case of an equator-on 
star as in our present case. We found a satisfactory fit of V and R light curves 
using two spots separated in longitude by about 140±30° and with radius of 
about r = 25±5° each, which corresponds to about 17% of the whole photo¬ 
sphere. Such a value refers to that component of spots unevenly distributed 
in longitude and, then, represents a lower value for the effective percentage 
of spotted surface. In Fig.[6j we plot the normalized flux in the B, V, and R 
filters with overplotted the fits from our spot modeling. The B model light 
curve is to some extent poorer with respect to the other fits, especially at the 
light minimum, where the observed light level is brighter than modeled. 


9 Conclusion 


Our investigation allowed us to measure the rotation period P = 4.43d of 
the dM component HD 155555C in the triple system HD155555. Comparing 
all components luminosities and effective temperatures with PMS models, we 
find that the hypothesis that all three components are coeval is sustainable 
with an age of 20±15 Myr, which is consistent with literature values and 
with the age of the f3 Pic Association. Then, the difference of rotation peri¬ 
ods can be explained considering the different nature of the components. The 
SB2 components that rotate with P = 1.687d have already reached a tidal 
synchronization, their orbital and rotation period being equal (Strassmeier et 
al. 2000; Cutispoto 1998). The single dM is still contracting and spinning up 
its rotation rate while approaching the zero-age main-sequence. The rotation 
periods of all components fit well into the period distribution of other /3 Pic 
members (see Fig. [7]). From the spot modeling we infer a high level of photo- 
spheric activity, with spots having a temperature contrast of 0.85 with respect 
to the unspotted photosphere and a covering fraction larger than about 17% 
of the whole photosphere. This level is consistent with the high activity level 
at the higher atmospheric levels reported in the literature. 
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